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Abstract. 

The projection formalism for calculating effective Hamiltonians and reso- 
nances is generalized to the nonlocal and/or nonhermitian case, so that it is 
applicable to the reduction of relativistic systems (Bethe-Salpeter equations), 
and to dissipative systems modeled by an optical potential. 

It is also shown how to recover all solutions of the time-independent Schro- 
dinger equation in terms of solutions of the effective Schrodinger equation 
in the reduced state space and a Schrodinger equation in a reference state 
space. 

For practical calculations, it is important that the resulting formulas can 
be used without computing any projection operators. This leads to a mod- 
ified coupled reaction channel/resonating group method framework for the 
calculation of multichannel scattering information. 
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1 Introduction 



In many applications, a quantum system of interest is part of a much bigger 
system, and the latter's state influences the system state. If the big system 
is represented by solutions of a Schrodinger equation in a big state space, it 
is desirable to find an effective Schrodinger equation in a small state space 
that describes how the small system of interest is affected by the embedding 
into the big system. Under certain conditions, a Schrodinger equation in the 
big state space can indeed be reduced to an effective Schrodinger equation in 
the small state space, in such a way that the interesting part of any solution 
of the full Schrodinger equation satisfies exactly the effective Schrodinger 
equation. 

Much effort has gone into solving this reduction problem, and in a sense 
it is well understood [1, 6, 13, 17, 28]. Exact expressions for the effective 
Hamiltonian can be given. In its exact form, the effective Hamiltonian is 
energy dependent and usually acquires a nonhermitian part; its eigenvalues 
describe the bound states and resonances of the reduced system. There are 
approximation schemes that compute the effective Hamiltonian (at least in 
principle) to arbitrary accuracy. 

Less known (but proved here) is that certain solutions of the full Schrodinger 
equation can be reconstructed from solutions of the effective Schrodinger 
equation, using little more than what is already available from the reduction 
process. 

The reduction is usually done for bound state calculations by the variational 
principle discussed in every textbook on quantum mchanics. For resonance 
calculations, the reduction may be done within the Feshbach [7, 9] projec- 
tion formalism (for an exposition see, e.g., KUKULIN et al. [17, Chapter 4]). 
In both cases, the reduced state space is finite-dimensional. For the calcu- 
lation of scattering states, the reduction may be done by means of coupled 
reaction channel equations, also called the resonating group method; a nice 
exposition is given in WiLDERMUTH & Tang [28]. In this case, the reduced 
state space is a direct sum of finitely many function spaces, one for each 
energetically admissible arrangement of particles into clusters, with states 
parameterized by coordinates or momenta of cluster centers only. 

The coupled reaction channel equations are numerically easy to handle, but 
the approximations involved in their derivation make an assessment of their 
accuracy difficult. On the other hand, the projection formalism gives in prin- 
ciple exact results, limited in accuracy only by the approximations made in 
the calculation of the effective Hamiltonian. However, this involves projec- 
tion operators, which are clumsy to use if an orthogonal basis is not easily 
available. 

Moreover, if the full Hamiltonian is already nonhermitian, or if the Hamilto- 
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nian is energy- dependent (such as in relativistic calculations; cf. [8, 14, 15, 
16, 12, 18, 19, 22, 26]), the standard derivation of the effective Hamiltonian 
is no longer valid. In the following, we shall remedy both defects. 

The paper is organized as follows. We first extend the projection formalism 
such that it apphes to the reduction of nonlocal systems and of dissipative 
systems modeled by an optical potential. Thus we work throughout with 
complex symmetric Hamiltonians, introduced in Section 2, and generalize in 
Sections 3 and 4 the traditional Feshbach projection approach to this more 
general situation. We then show in Section 5 that the formalism in fact 
allows to recover all solutions of the time-independent Schrodinger equation 
in terms of solutions of two Schrodinger equations, one in the reduced state 
space and the other in a reference state space. The resulting formulas are 
closely related to those of time-independent perturbation theory. 

For practical calculations, the resulting formalism is revised in Section 6 so 
that it can be used without computing any projection operators. Section 7 
shows how the use of doorway operators gives flexible approximation schemes 
for the exact formulas derived in Section 6. A backward error analysis dis- 
cussed in Section 8 allows the estimation of the reliability of approximate 
solutions of Schrodinger equations obtained by this or any other method. 

Specific choices of the embedding map lead in Section 9 to a modified cou- 
pled reaction channel/resonating group method framework for the calcula- 
tion of multichannel scattering information, in which solutions of the full 
Schrodinger equation can be obtained from solutions of coupled reaction 
channel equations for the effective Hamiltonian. In principle it is capable 
of arbitrarily accurate approximations to the full dynamics, and shares this 
feature with the two Hilbert space method of Chandler & Gibson [3, 4], 
which partly inspired the present investigations. 

The theory is presented in a fully rigorous manner, allowing for unbounded 
operators by using in place of Hilbert spaces a pair of dual topological vector 
spaces. The most useful results are in Sections 6, 7, and 9. 

The reader interested in the results but not in mathematical rigor may 
omit all references to spaces and topology, may think of all spaces as finite- 
dimensional and of operators as matrices, and may skip all proofs. In particu- 
lar, of Section 2 introducing the basic terminology, only (l)-(3) are essential, 
and most of Sections 3-5 can be skimmed. 



2 Symmetric operators 

Lin(V, W) denotes the space of continuous linear mappings between two topo- 
logical vector spaces V and W, V* = Lin(V, C) denotes the dual space of 



3 



continuous, complex- valued linear functional on V, and LinV = Lin(V, V) 
denotes the algebra of continuous linear transformations of V. 

In the following, HI is a complex topological vector space with a definite, 

continuous symmetric bilinear form, providing a natural embedding of H 
into the dual space H*. We refer to H* as a state space, and to the ip eM* 
as states. We write the pairing and the bilinear inner product as 

(p'^^P = ^'^(p for e H, -0 e H*. 

The notation is chosen such that it looks as closely as possible like standard 
finite-dimensional linear algebra. 

We say that a sequence (or net) ipi e H* converges weakly to -0 e H*, and 
write ipi ^ ip, ii 

lim (f^ipi — (p^il>\/(p e H. 

l—fOO 

We extend the bilinear inner product to arbitrary pairs {(p, ^) e H* x M* for 
which 

ip^i/j — lim (fjipi 

l—fOO 

is defined and independent of the weakly converging sequences (or nets) 
pi (p, ipi tp of (pi,tpi G H. (In the applications, this allows to form the 
inner product of state vectors corresponding to bound states and resonances 
but not that of scattering states.) 

Complex conjugation is denoted by a bar, and, with the notation ip* — 'ijj'^, 
the Hermitian inner product on H is 

{(f\ip) = (f*ip for (p,ip 

The associated Euclidean norm is 

IIV'II = Vi^, 

and Si is the closure of H in H* with respect to the Euclidean norm. Thus 
H C H C H*, and H is a Hilbert space. 

In the applications, HI is a space of sufficiently nice functions (namely arbi- 
trarily often differentiable, with compact support) on some finite- or infinite- 
dimensional manifold, the inner product of two functions is some integral of 
their pointwise product induced by a nonnegativc measure on the manifold, 
and H C H C H* is a Gelfand triple (or rigged Hilbert space). (For these 
concepts, see Gelfand & Vilenkin [10], Maurin [21]. For an exposition 
in physicists' terms see KUKULIN [17, Appendix A]; cf. also BOHM [2].) 

The transpose of a linear operator A G Lin(lHI, H*) is the linear operator 
e Lin(H, H*) defined by 

{A^ififi/j = ip^Ail; for -0 G H, 
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A* — defines the adjoint of A, and 

ReA=-iA + A*), ImA= -(A- A*) 

define tlie real and imaginary part of A. Clearly, {ABY = B^A^ and 
(AB)* = B*A*. The operator A G Lin HI is called symmetric if A^ — A on 
H, Hermitian ii A* = A on H, and positive semidefinite if 

ip*Aip > for aU ip eU. 

In particular, A is Hermitian if and only if ImA = 0. We extend symmetric 
operators A e Lin HI to Lin H* by defining e H* for e H* by 

if^A^} = {Aip)'^^y(p e H. 

In the following, L e Lin HI is always a symmetric operator such that 

Im L is positive semidefinite. (1) 

(In particular, this includes the case where L is Hermitian since then Im L = 
0.) 

Since i/j*Li/j = ip* (Re 1)1/; + iip* {Im L)tl; and both V*(ReL)V^ and il;*{ImL)tp 
are real, (1) is equivalent to 

lmi;*Lijj>0 forall^eH. (2) 

In particular, the spectrum of L is in the complex upper half plane. If L has 
a spectral resolution then, since |A + ie] > £ for Im A > 0, we conclude that 
{L + i£)~^ exists for all £ > as a bounded operator on EI with spectral norm 

IKL + ie)-^ < foralle>0. 

The traditional situation is the one where L = E — H with a complex energy 
E satisfying Im E' > and a Hamiltonian H — Hg — ^V . Here Hs,T e Lin H 
are symmetric and Hermitian, and F is positive semidefinite. In the most 
important case of a conservative system, F = and H is Hermitian. However, 
care is taken that all our results hold in the nonhermitian case, corresponding 
to dissipative systems with an optical potential that contributes to F. 

The use of L helps to avoid a multitude of expressions involving E — H 01 
E — H + ie. Since the (time-independent) Schrodinger equation H'ip = Eip 
takes the simple form 

= 0, (3) 

this makes the formal manipulations independent of energy and free of ref- 
erences to the Hamiltonian H, and thus much more readable. (Of course, in 
actual calculations, E and H reappear.) 

More generally, (3) also covers nonlocal problems with a nonlinear depen- 
dence of L on E, and therefore can be used for the Bethe-Salpeter equations 
arising in bound state and resonance calculations for relativistic systems 
[14, 15, 18, 19, 26]. 
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3 Reduction of the state space 



Let H and Hcff be topological vector spaces with a definite, continuous bilin- 
ear inner product, related by the embedding map P, an injective, closed 
linear operator from E[*g to H* satisfying P = P and 

Pipes e H for aU IpeS e Heff. 

We want to relate a Schrodinger equation Lip = in the full state space H* 
to an effective Schrodinger equation LesipeS = in the reduced state space 

P* = maps B* to H*g and H to Heff. Moreover, P*P : H*g H*g is 
invertible since P is closed and injective. The pseudo inverse 

pi ^p*pyip* 

maps H* to H*g and possesses the properties 

{PP^f = PP^, P^PP^ = P^, P^P = 1. (5) 
This imphes that 

g := 1 -PP^ = 1 - (P^)^P^ e Lin H* (6) 

satisfies 

P'Q = P^Q = 0, QP = 0, Q^ = Q^ = Q. (7) 

and hence is the orthogonal projection to the orthogonal complement of the 
range of P. 

A Q-resolvent of a symmetric operator L e Lin HI is a symmetric operator 
G e Lin HI satisfying 

GP = 0, GLQ = Q. (8) 
Since GQ = G{1 — PP^) = G, the symmetry of G implies 

QG^G^ GQ, QLG = Q = GLQ. (9) 

Formally, G — Q{QLQ)~^Q, but {QLQ)~^ is only defined on the range QW 
of Q. 

The following discussion generalizes the Feshbach projection formalism which 
is obtained in the special case where L — E — H and H is Hermitian and G is 
the ordinary resolvent of QHQ in QW (cf. the development in KUKULIN et 
al. [17, Chapter 4]; in their notation, G — Gq{z) is called the 'orthogonalized 
resolvent'). 

3.1 Proposition. 

(i) If L has a Q-resolvent G, it is uniquely determined, the operator 

Les P'^LP - P'^LGLP = P'^{L - LGL)P = P^ L{P - GLP) (10) 
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is symmetric, and Im L^s is positive semidefinite. 
(a) For arbitrary ip^s G E[*g, the vector 

(P-GLP)V^effeH* (11) 

satisfies 

P^LV; = LeffV'eff, QLi^^O, Veff^^V- (12) 

Proof. If (9) holds with G' in place of G then 

G" = G"Q = CQLC = G'LG = CLQC ^QG^G, 

giving uniqueness. The vector (11) satisfies 

QL^ = QL{P - GLP)^^s = (<5 - QLG)LPi;^s = 0, 

P^LV- = P^L(P - G'LP)^eff = L^si^es, 
and, since P^G = P^QG = 0, 

This proves (ii). Since P — P,we have 

= ^*(P^)*Leif^eff = ^eff-^eff^eff- 

Thus, for arbitrary ■0efr £ Heff' Im V'^-^'effV'eff = lm.'4>*L'^ > 0. This 

proves (i). □ 

The second term 

A := P^LGLP (13) 

in the definition of L^s is called the optical potential induced by the re- 
duction process. (The name is explained in Taylor [27, p. 385].) 

In the special case where P^P = 1 and L = E — H, G — G{E) and hence the 
optical potential A (P) = P^{E—H)G{E){E—H)P is energy-dependent (and 
nonlocal) and we have L^s = E — H^q[E) with the effective Hamiltonian 

H,^{E)^P'^HP + /:^{E). 

Thus the optical potential causes energy shifts in the eigenvalues of the pro- 
jected Hamiltonian P^HP. We also note that the reduced Schrodinger 
equation is generally a nonlinear eigenvalue problem 

H,^{E)ij = P^. 

Frequently, the energy-dependence is ignored; however, nonlinear eigenvalue 
problems for nonlocal Schrodinger equations arising from Bethe-Salpeter 
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equations for relativistic problems were actually solved, e.g., [14, 15, 18, 19, 
26]. 

Note that the resonating group method (Wildermuth & Tang [28]) works 
with coupled reaction channel equations derived from a projected Hamilto- 
nian and hence misses the optical potential; an energy-independent term is 
added instead on a phenomenological basis [28, Chapter 8.2]. An alternative 
exact method is derived below. 

We return to the general case. 

3.2 Proposition. The following identities hold: 

L,sP' = P^(l - LG)L, {PYL,s = L{P - GLP), (14) 
and, if L and L^s are invertible, 

P'^L-^P = P^'PL-^P'^P. (15) 

Proof. Since (L - LGL)Q = L{Q - GLQ) = 0, we have 

LeffP^ = P^(L - LGL)PP^ = P^(L - LGL){1 - Q) ^ P'^{L - LGL). 

This gives the first equation in (14), and the transpose gives the second 
equation. For invertible L, Leff, we conclude from (14) and GP — GQP — 
that 

L^sP^L-^P = P^(l - LG)P = P^P, 
hence P'^PLJP'^P = P'^PP^L~^P = P'^L-^P. □ 

(15) implies that for invertible L and Leg, the part of accessible from H*^ 
can be computed from the knowledge of Lgfj alone. Similarly, our next result 
says that, if L has a Q-resolvent , all solutions of the Schrodinger equation 
Lil) = can be computed from solutions of the reduced Schrodinger equation 
-^eff^eff = and a knowledge of the correction operator 

R^GLP (16) 

occurring in (10) and (11). 

3.3 Theorem. Suppose that L has a Q-resolvent G. 

(i) For any tpes £ with LesipeS — 0, the vector tp e M* defined by (11) 
satisfies Lip — and 

V'eff^PVeH.V (17) 

(a) For any ijj G H* with Lip = 0, the vector (17) satisfies Legtpes — 0, and 
we can reconstruct ip from (11). 
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Proof, (i) Multiplication of the second equation of (14) with V^eff gives 

= (P^^LeffTAeff = L{P - GLP)^,s = Lij, 
and (17) follows from Proposition 3.1 (ii). 
(ii) Multiplication of the first equation of (14) with ijj gives 
Le#eff = L,sP'^ = P^(l - LG)Lij = 0. 
Since P^^eff = PP^i^ = (1 - Q)i^ ^i^-Qi^, (11) follows from 

{l-GL)ij = {1 - GL)ij - {Q - GLQ)i; = {1 - GL){1 - Q)i; 
= (1 - GL)P^P,s = (P - GLP)i;,{i 

since 

GLtjj = 0. (18) 

n 

Note that, while LcsipcS = looks like a Schrodinger equation, this is gen- 
erally a nonlinear eigenvalue problem. Indeed, if L = L[E) = E — H then 
Les{E) is generally a nonlinear analytic function of E, and the effective 
Hamiltonian 

Hes{E) :— E — Les{E) 
has a nonlinear dependence on the energy. 

In many cases of interest, the Q- resolvent does not exist for the Hamiltonian 
H of interest. But frequently the Q-resolvents G^ for the perturbed Hamil- 
tonians H — is corresponding to — L + ie exist for all e > 0, and are given 

by 

G, = Q{QLQ + te)-^Q. 

If the limit 

P :=limG^ L^P (19) 

exists, most of the preceding proof still goes through, with R in place of GLP 
and in place of P^LG, and the effective Hamiltonian (usually) acquires 
a nonhermitian part. 

The only exception is the second half of statement (ii), which must be mod- 
ified. (This is most conspicuously seen when P = 0, where L^s = and 
the reduced Schrodinger equation provides no information at all.) Inspection 
of the proof shows that (18) fails. Thus, if ■0 e H* satisfies Lip — 0, the 
expression 

i;^:=tP-{P- R)^,s (20) 

need not vanish (as predicted by (11) under the stronger assumptions), but 
only the much weaker equation P'^i/j-^ — follows. Thus the reduction does 
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no longer allow one to recover all solutions of Li/j = 0. By Theorem 3.3(i), 
t/j-^ is also a solution of the Schrodinger equation in H*, and since P'^il) — 
implies ipgs = and hence tp-^ = ip we have 

{^^ \ LiP ^ 0} ^ {iP \ LiP ^ 0, P'^i) = 0}. 

We discuss later (Theorem 5.2) how to access this missing part. 



4 Perturbation theory 

The computation of the Q-resolvent is traditionally done using perturbation 
theory. It is assumed that the inhomogeneous Schrodinger equation for a 
related reference problem is explicitly solvable, and the problem of interest 
is considered as a perturbation of the reference Schrodinger equation. 

We therefore assume that we have a symmetric operator Lref G LinM satis- 
fying 

LreiQ = QLrch (21) 

this commutation relation can always be achieved by taking an arbitrary 
symmetric approximation Lq to L and putting 

L,ef = Lo - (1 - Q)L^Q - QLo{l - Q). (22) 

Let G be a Q-resolvent of L and put 

V = Lref -L, T = V + VGV, n:=l + GV. (23) 

V is called the interaction, T the (Q- version of the) transition operator or 
T-matrix, and Q the (Q-version of the) Moller operator. This terminology 
is justified by the close formal relations of the properties of T and fl, derived 
below, with those of the T-matrix and the Moller operators of standard 
scattering theory. Indeed, for P = 0, Q = 1 (which is uninteresting from the 
point of view of effective Schrodinger equations), the results reduce to those 
of standard perturbation theory. 

4.1 Proposition. If (21) holds then 

nP={l- GL)P, n^P = P, (24) 

T^Vn^ fFV, (25) 

P'^mP = Leff = P^(Lref - T)P. (26) 

Proof. Since 

GLrefP = GQL^^iP = GU^iQP = 

we have 

GLP = GL,^iP - GVP = -GVP, 
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P'^LG = {GLPf = -{GVPf = -P'^VG, 

hence 

= P + Gyp = (1 - GL)P. 

Together with n'^ P = (1 + VG)P = (1 + yG'Q)P = P, this gives (24). (25) 
follows directly from (23). Since 

P^LQP = P'^L{1-GL)P^ L^s 

= P^LP - P^LGLP = P'^LP + P^LGVP 

= P^(L,ef - y)P - P^VGVP = P^(L,ef - T)P, 

(26) follows. n 



4.2 Proposition. If fl is invertible then 

Gref = ^^-'G (27) 

is a Q-resolvent of Lref, and 

Q = (1 - M^)-\ wiere = G^efV. (28) 

Proo/. 

0(1 - W^) = - f^GrefV ^n-GV ^1, 

gives (28). Since G^efP = Q'^GP = and 

GrefLr^fQ = GrefLQ + Gref^Q = Q'^GLQ + W^Q = (1 - H/)g + Vl/Q = g, 
Gref is a g-resolvent of Lref. O 

4.3 Theorem. Suppose that Gj-cf is a Q-resolvent of a symmetric operator 
Lj-cf satisfying (21), and suppose that Q with (28) exists. Then G — OGref -js 
a Q-resolvent of L = L^ef — y, and with 

T^V + VGV, (29) 

(23)-(27) hold. Moreover, we have 

n^i + wn^i + nw, (so) 

T = V + TW = V + W^T, (31) 
G = Gref + H/G = Gref + GM^^. (32) 

Proof. (30) follows directly from fi(l — W) = (1 — W)Q = 1 and implies 
Q = 1 + QGref^ = 1 + GF, hence (23). Since GP - QG^efP = and 

GLQ = GL,^iQ - GVQ = nG,^iL,^fQ - GVQ = QQ - GVQ = Q, 



11 



G is a Q- resolvent of L. Thus Proposition 4.1 applies, and gives (24)-(26). 
(27) is obvious. The first equality in (31) follows from (25) and (30), and the 
second follows by transposing the first equation. The first equation in (32) 
follows from 

WG = G^eiVG = G^eiQVQG = Gref(<9^ref(5 - QLQ)G = G- G^ef, 

and the second follows again by transposing the first equation. □ 

Inserting the definition (28) of W into (32) gives the (Q-version of the) 
Dyson equation 

G = Gref + Gj-eiVG = Gref + GV^Gj-gi- (33) 

If the spectral norm oiW — G,-e{V is smaller than one, (30)-(32) can be used 
to calculate iteratively the MoUer operator Q, the transition operator T and 
the Q-resolvent G. To lowest order, we get 

Q « 1 + GrefV, T^V, G^ Gref + G,ef1^G,ef • 

When Heff = {0}, this is the Born approximation, and when is finite- 
dimensional, this is a version of the distorted wave Born approximation 
(see, e.g., Newton [24, Section 9.1]). In the Born approximation, we simply 
get 

in second order, 

Leff ^ P'^LP - P'^VG.efVP, (34) 

giving the approximation A ^ P'^VG^eiVP for the optical potential (13). 
Further iteration gives the Born series 

T^V + VG^efV + VG^eiVG^efV + ..., 

giving the exact optical potential 

A = P^yGrefV^P + P^VG,e{VG,e{VP +.... 

Thus we have recovered a generalized version of traditional perturbation 
theory. For M^s = {0}, P = 0, Q = 1, the equations obtained above reduce to 
those of standard perturbative scattering theory; the missing part - that one 
gets the scattering solutions of the Schrodinger equations - follows in the next 
section. Of course, from the point of view of effective Schrodinger equations, 
the case P = is completely uninteresting; however, it is instructive in that 
it shows that the methods used to solve scattering problems apply with small 
modifications to the problem of finding effective Hamiltonians. 

Similarly, if Heff is the eigenspace of if^ef corresponding to a bound state of 
the reference system and P the orthogonal projector to this space, we get 
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the situation leading (in the Born approximation) to Fermi's Golden Rule; 
cf. KUKULIN et al. [17, Section 4.4]. For a nondegenerate bound state, Hefj 
is one-dimensional, and again the point of view of effective Schrodinger equa- 
tions is empty. For degenerate bound states, however, we recover a nontrivial 
low-dimensional eigenvalue problem as effective Schrodinger equation. 



5 Solving the Schrodinger equation 

We are now ready to express all solutions of the Schrodinger equation in terms 
of the solutions of a reduced Schrodinger equation and special solutions of a 
reference Schrodinger equation. 

5.1 Proposition. If 

= PV, ^ref = Q(l - W)ilj (35) 

then 

P^^ref =0, V = ^^(^'V'efr + Q^ref) • (36) 

If (36) holds then 

Weff = P'^O^LV', (37) 

^refV'ref = QL'^, (38) 
Lll; = (P^)^LeffV'eff + Wref " (1 " Q)^'' L,^^,^^. (39) 

Proof. Since 

QW = QGrefV = G,,iV = W, 

(35) implies 

= (1 - PP^ - W)i^ = (1 - W)i^ - Pipes, 

hence 

^{Pi^oS + Q^Prci) = ^(l - W)iP = V, 
P'^^.ef = P^Q{1 - W)ij = 0. 

This gives (36). Now suppose that (36) holds. By (24), 

P'^n^L = {nPfL = P^(l - LG)L = P^L(1 - GL), (40) 

and 

P'^n^LQ = P^L(Q - GLQ) = 0. (41) 
Since (1 - QW)n = {1-W)n = l, this implies 

P'^n^m = P^Q^L(l - QW)n = P^Q^L. (42) 
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Using (36), (42), (41), and (40), we find 

= P'^n^LPlP.s = P^i^(l - G'L)PV^eff = Weff, 

giving (37). Since QLj-f^^W = QLrefGref^ = QV, we have 
hence 

QLV' = QLVL{PiI)^q + Q^/^ref) = L^eiQiPtpeS + QV'ref) 

by (35), giving (38). Finally, X := 1 - (1 - Q)Q^ satisfies 

XP = P - (1 - Q)O^P = P - (1 - g)P = QP = 
by (24), hence 

XQL = X{1 - PP')L = XL = (1 - (1 - g)Q^)L = L - {P'fP'^^FL, 
so that 

LV' = ((P^^P^n^ + Xg)LV; = {P'fL.si^.S + ^ Wref 

by (37) and (38). □ 

Since the formulas in Proposition 5.1 and L^e = P^LflP do not involve G, 
we can take limits and obtain: 

5.2 Theorem. Let (jref,£ be a Q-resolvent of I/ref,e with limLref,£ = Lref- 
Suppose that 

V^^ := limGref,£^ and n = {1 - W)'^ 

exist. If 
then 

P^^ref =0, V = ^^(i'^eff + Q^ref) • 

Moreover, with L^s :— P^LQP, we have 

6 State space reduction without projections 

Projection operators and the associated Q-resolvents are clumsy to use if 
an orthogonal basis is not easily available. We therefore revise the above 
formalism so that it can be used without computing any projection operators 
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or Q-resol vents. Additional flexibility is gained by using in place of the 
reference operator Lj-^f G Lin EI an appropriate operator Pq : Mq ^ H* from 
the dual of a reference space Mq. For exact solutions, this reference space 
must be at least as big as H; however, in Section 7, we choose Hq to be 
a smaller space in which numerical calculations are tractable, and obtain 
practical approximation schemes for the correction operator and the effective 
Hamiltonian. (This is related to the two Hilbert space method of CHANDLER 
& Gibson [3, 4], and indeed was inspired by their work.) 

As before, L G LinH is assumed to be symmetric, and P : IHI*^ H* is 
assumed to be a closed injective linear mapping satisfying P — P and 

Pipes G H for aU -063 G Heg. 

Let Pq -.Mq ^ H* be a closed, surjective linear mapping satisfying 

PoV^o e H for all V'o e Hq. 

Then P^ maps H* to and M to Hq. Moreover, PqP^ : M* ^ H* is 
invertible since Pq is closed and surjective. 

The pseudo inverse P^ of the injective P satisfies as before 

p^ = {p'^py^p'^^ p^p = 1, p^pp^ = p^ 

but the pseudo inverse Pq of the surjective Pq satisfies 

Pq^Pq{PoPo) \ -Po-Pq = 1, PqPgPq ^ Pq ■ 

6.1 Proposition. If there are linear mappings R : E[*g — > H* and L' : E[*g — > 
sucii tiat 

PVeff e M for all ipes e IHeif ; -L'-^o e Meg for all ipo £ Eleg 

and 

Po^LP + P^PL' = Pq^LP, P^P = 0, (43) 

tiien 

Leff P^LP - P^LP G Lin H (44) 
is a symmetric operator satisfying 

L,s = P'^LiP - R), {PYL,s = L{P - R). (45) 

Moreover, if L and L^s are invertible then 

P'^L-^P = P'^PLjP'^P. (46) 

Proof. Symmetry is obvious. Multiplication of (43) with Q{P^y gives 

QLP = QLR + QPL' = QLR. 
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Now 

R'^LR = R'^L{Q+{PYP'^)R = R^LQR 

= {QLRfR = {QLPfR = P'^LQR = P'^LR 

since QR={1- {P^YP'^)R = R, hence 

Leff = P'^LP - R^LR = P'^LP - P'^LR = P^L(P - R). 

This is the first half of (45). By multiplication with (-P^)^, we find 

(PT^eff = {PTP^L{P-R) = {l-Q)L{P-R) 
= L{P -R)- QLP + QLR = L{P - R) 

giving the second half of (45). If L and Les are invertible then (45) implies 
hence 

P'^PL-^P^P = P^L-\P^fP'^P = P^L-^P, 
giving (46). □ 

It is not difficult to see that R and L^s from Section 3 are an instance of this 
construction, with EIq = EI and Pq = 1- 

6.2 Proposition. 

(i) For arbitrary if) eM* , the vectors 

V^eff = e h:^, Vo = n'(V'-(^'-i?)V'eff)eH* (47) 

satisfy 

P^Po^o = 0, V = - R)i^e{i + Po^o. (48) 
(a) If (48) holds for some tpes e lHI*g, e Hq then 

Les^es =iP- RfLtP, (49) 

P^LPoi/jo = {Po -{P- R)P'PofLiP, (50) 
Lil; = (P^)^Leff^eff + (P^fiP^LPoi^o). (51) 

Proof, (i) follows from 

PoV'O = PO^O (V' - - R)M = ^ - (P - P)Veff, 
P^Po^O = P^ij-P^{P-R)llJen=P^i^-P^Pi^eS 

= P^i/j - P'^PP^ip = 0. 

(ii) Since by (43) 

P^LV; = R^LPoiJo = {P;fLRfijo = {P;fLP~P^PLY^o 
= P^LPo^o - L'^P^Poi^o = P^LPoiJo, 
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(49) follows from (45) and (48): 

= P^LV; - P^LPoV'o = P'^Li; - L'^ = (P - RfLi^. 

(50) holds since by (48), (45) and (49), 

Po^LPo^o = P^Lij-Pi'L{P-R)'i/;,s = PifLij-P;[{PYLesAs 
= P^LiP - P^{PY{P - RfLiP 
= {Po-{P-R)P'PofL,l;. 

(51) holds since by (47) and (45), 

{PlfP^LP^il^o = LPoiJo = LPoP^{iJ^iP-R)iJes) 

^ _ (p _ p)^^g) =Li;- (P^)^LeffV'efr- 



We now have the following projector- free and Q-resol vent-free version of The- 
orem 5.2, constructing all solutions of the Schrodinger equation Lip = in 
terms of two simpler Schrodinger equations. 

6.3 Theorem. 

(i) If Ira L is positive semideGnite then Im L^s is positive semideGnite. If also 
Pq = Pq then also Im P^LPq is positive semidefinite. 

(ii) For arbitrary V' e M*, V'eff e 11*^, ipo e satisfying (47) or (48), 

In particular, to find all solutions of Li/j = 0, it sufEces to solve the two 
problems 

P^LPoi^o = 0, P^PoV'o = 0. 

Proof, (i) For arbitrary ipcS ^ Hcffi define ip by (48) with ipQ = 0. Then 
P^tp = P^{P - R)'ip^s = P^Pip^s, hence tp^s = P^ip- Now (51) gives 

Li^ = (P^^ Weff, 

and since (P^)* = {P^Y , we get 

riip = riP'YLoS'ipcS = {P''il>yL,siJ,s = CfrWcfr. 

Hence ip*ff{lm Lcs)'ipcH = ip*{lTaL)ip > 0, and ImLcg is positive semidefinite. 
If also Pq = Pq then Pq* = Pq and Im PqLPq is positive semidefinite since 

^*(ImPo^LPo)^o = (Po^o)*(ImL)(Po^o) > 0. 

(ii) The forward implication follows directly from (49) and (50), the reverse 
implication from (51). □ 
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6.4 Remark. If = then <^^V = vlsiP'^P + R^R)As, so that 

GeS ~ P + FlF R 

is the effective metric induced on Hcff. Note that it is generally not the 
original metric in Heif, not even when P is an orthogonal projection. 



7 Form factors 



In this section we show how to obtain efficiently approximations to the cor- 
rection operator R. This leads in Section 9 to a modified coupled reaction 
channel/resonating group method framework for the calculation of multi- 
channel scattering information. 

We emphasize that the formulas derived in this section no longer involve a 
pseudo inverse. In particular, they can be used even when P is not injective 
and Pq is not surjcctive. (Hovever, since the assumptions under which the 
formulas arc derived are then violated, they lose some information and hence 
give only approximate effective Hamiltonians.) 

7.1 Theorem. Let 

Lo:=P^LPo, P,:^P'^Po, U := P^ LP - P^ PL (52) 
with a symmetric L e Lin Hes (and hence in Lin IH[*gj, and write 

Ls PI 



. , e Lin (Ho © Heff) with L^^ Lo + ie. 
Pi —le ' 



(i) If the strong hmit 



exists then (43) is solved by 

R^PoFo, L'^L + Fi. (54) 
(a) Relation (53) holds with 

Fo = lim Foe , Fi = hm Fi, , (55) 

where 

Fu = (PiL:^Pi + ie)-'PiL;% Fo, = L;\U - P^Fu), (56) 
if these limits exist. 
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Proof, (i) follows from 



P^LR + P^PL'\ ^ fP^LPoFo + PJPF, + PjPL 
P^R J V P^PoFo 



Fo\ , fP^PL 



FJ V 



by looking at the upper and the lower part separately, (ii) follows from the 
equation 

which is easily verified by substitution. □ 



In principle, L in (52) may be chosen arbitrarily. However, for numerical 
calculations it may be advisable to choose L in such a way that U (which 
replaces the interaction V in the projection approach) becomes small in some 
sense. This has the beneficial consequence that then the numerical approxi- 
mation errors have a much smaller effect on the calculated solution. 



In practice it is impossible to compute the exact correction operator and 
hence the exact effective Hamiltonian, since these tend to be exceedingly 
complicated. One therefore exploits physical intuition to select a space Hq of 
manageable complexity whose dual contains the doorway states believed 
to mediate the interaction of the reduced system and the unmodelled envi- 
ronment. Hq is embedded into M* by means of a doorway operator Pq 
that is now no longer surjective. Fortunately, the formulas (43) and (44) 
defining R and L^s do not depend on pseudo inverses, and hence make also 
sense in this case. (43) and (44) now only yield approximate solutions for 
the correction operator and an approximate effective Hamiltonian. However, 
these approximations become better and better as the range of Pq covers a 
bigger and bigger part of H. 

The situation is fully analogous to numerical discretization schemes that are 
necessary to solve all but the simplest partial differential equations; the only 
difference is that in the present context it frequently makes sense to consider 
approximations in manageable function spaces, so that one does not discretize 
completely. 

A proper choice of the doorway operator Pq makes the computations more 
tractable. At the same time, it limits the formal complexity of the optical 
potential 

A = R^LR = F^{P^LPo)Fo, (57) 
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the second term in Lgff. As one can see, Pq specifies the allowed form of the 
optical potential while the form factor Fq specifies the coefficients in the 
optical potential, and thus introduces energy- dependent running coupling 
constants. 

The art in applying the reduction technique consists in finding embeddings P 
that 'dress' the subsystem of interest in a sufficiently accurate way, a doorway 
operator Pq embedding the relevant doorway states, and an operator L such 
that U is small, and the limit (53) exists and can be approximated efficiently. 

Setting Pq = gives = 0, A = 0, and hence the trivial approximation 

Leif !^ P^LP. (58) 

This simply amounts to discarding the interaction of the subsystem with 
the rest of the system. The choice Mq = M, Pq = P is not better since 
then P^PFq = P'^PqFq = P^R = by (43), and hence Fq = 0, i? = 0. 
This is not surprising since we expect that the doorway operator Pq should 
incorporate additional information about doorway states not yet represented 
in the subsystem but significantly interacting with it. 

If Ho and HIeff are finite-dimensional then (52) defines finite-dimensional ma- 
trices, and the computation of the form factor amounts to solving the matrix 
equation 



p. ; w vo ' '''' 

with a complex symmetric (and for L = L*, Pq = Pq Hermitian) coefficient 
matrix. (59) can be solved efficiently by sparse matrix methods (cf. Duff 
et al. [25, 5]) if suitable localized basis functions are used to construct P and 
Po. 

In practice, it may be useful to employ in combination with discretization 
methods a complex absorbing potential in place of the +ie. In particular, if 
one proceeds as in Neumaier & Mandelshtam [23] one gets a quadratic 
eigenvalue problem that can handle all energies in a certain range simulta- 
neously using harmonic inversion (MANDELSHTAM & TAYLOR [20]). 

If Hcff is finite-dimensional but Hq is a function space then Lq is a differen- 
tial or integral operator on Mq. By solving suitable differential or integral 
equations we can find the vector- valued functions 

Bq := hm LZ^P^LP, Bi := lim L'^PIP. 

and the complex symmetric matrix 

Gi := Pi El = lim(PiL7^Pi^ + ie). 

Noting that L' — Fi'ii L — 0, the formula (55) for the form factor becomes 
Fq = Bq- BiL', where L' = G^^PiBq. 
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8 The quality of approximate state vectors 



In practice, it is usually impossible to find exact solutions of a Schrodinger 
equation Ltp = 0. On the other hand, in real applications, L is never precisely 
known either. Hence it makes sense to assess the quality of an approximate 
state vector ijj by trying to modify L a little to an operator L that satisfies 
Lip — exactly. If the modification L — Lis within the accuracy to which L 
is known, we are confident that is a good approximation to the true but 
unknown L. 

This way of assessing the quality of an approximate solution of a problem is 
widely used in numerical analysis (see, e.g., Wilkinson [29]) and is known 
under the name of backward error analysis. Here we give a backward 
error analysis for the equation Ltp = 0, and deduce guidelines for quality 
assessment of approximate state vectors. 

8.1 Theorem. Let L G Lin HI he symmetric. Then, for arbitrary ip & M* 
with finite ip'^ip ^ 0, the modified operator 



satisfies Li/j = and 

tT{L-Lf^TM, (61) 

where 

Proof. Since the formulas are invariant under scaling ip we may assume that 
ip is normalized to norm 1. Then 

ip'^ijj = 1, V^LV' =: A, V^L^V =: /x, 
with real A, The operator 

A := L - L = iptp'^L + LipilF - XipijF 
satisfies A'^ = Lip, hence Lip — 0. Since 

/\LiP = XLiP \^)ip, 

we find 

A^ = LiPip'^L +{fx- X^)i^iF. 

A^ maps H to the two-dimensional space spanned by and L'^, hence is 
trace class. Using the formula tr if^"^ — ip'^(p, we find 

tr A^ = /i + (yu - A^) = 2// - A^ = tl{iP), 
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giving (61). n 

In the conservative case where ImL = 0, and hence L is Hermitian, it is 

possible to show that the choice (60) is best possible. Note that a state 
vector ip &W with ijj = ip and finite i/j'^ip is now in the Hilbert space H. 

8.2 Theorem. Let L e LinH be Hermitian, and suppose that ip e EI\{0} 
satisfies ■^^ — and Lip e H. 

(i) Any symmetric and Hermitian L e Lin H with Lip — satisfies 

ti:{L-Lf>TL{ip). (63) 
EquaUty in (63) is achieved precisely when L — L. 
(a) We always have 

(LiPfLib (LiPfLiP , , 

Jfl particular, tl{iP) = if aud oiiiy if L-i/' = 0. 

Proof. By the preceding theorem, L = L gives equality in (63), and is a good 
choice since Lip — 0. Hence suppose that L ^ L. Without loss of generality, 
(L — LY is trace class (otherwise the trace is infinity and (63) is trivially 
satisfied). Since (L — L)'0 = Lip — Lip = we have 

tr(L-L)(Z-L) =trA(Z-L) 

= iP'^L(L - L)iP + iP'^(L - L)LiP - XiP'^(L - L)i) = 0. 

Therefore 

tr(L - Lf - TLiip) = tr(Z - L)2 - tr(L - L)2 
= tr(L + L-2L)(L- L) 
= tr (L - L)2 - 2 tr(L -L){L- L) 
= tr(Z - L)2 > 

since L ^ L. Therefore, (63) holds ior L ^ L with strict inequality. This 
proves (i). Since ip — ip, the second term in (61) is nonnegative. This gives 
the upper bound in (64) and implies the final assertion. The lower bound 
follows from the Cauchy-Schwarz inequality [ip'^Lipy < ip'^ip- [Lip)'^Lip. □ 



(60) is the orthogonal projection of L to the orthogonal complement of ip, 
and a/ tl('0) measures, in a sense, its deviation from L. Therefore, Tiiip) (or 
its square root) serves as a useful measure for the quality of an approximate 
solution ip of the Schrodinger equation with Hermitian L. 
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In the nonhermitian case, it seems possible that ti{iIj) = even if Lt/j ^ 0. 
(A finite-dimensional example is L = g), -0 = (J) which has = and 
= for a = ^/2, and for a = 1 even ri{ip) < 0.) Thus, unless L* = L, 
the measure tl{iP) might be sometimes too optimistic. However, in practice 
L is nearly Hermitian and the use of tl{iI^) should cause no problems. 



9 Multichannel scattering 

We now apply the above to the multichannel approach discussed below. A 
projection approach to multichannel scattering leading to effective Hamil- 
tonians is discussed, e.g., in Newton [24, Section 16.6], but the equations 
derived there appear not to be suitable to numerical approximation. A more 
useful formulation is given by the present equations from Theorem 7.1, with 
P and Pq as given below. 

An arrangement is a partition A of the system of particles into clusters 
i & A, with correct assignment of distinguishability. At a fixed energy E, 
those arrangements arc relevant that contain channels defined by cluster 
bound states with energies Ej, such that 

^Ei<E + AE 

i€A 

where AE is zero or a small quantity. These open or nearly open channels 
are assumed to correspond approximately to states in a n^-dimensional space 

with basis functions 

(pAk{x) =Y[^iii^{xi), (65) 

ieA 

where Xi is the vector of coordinates of particles in cluster i, and the (pii{xi) 
are translation invariant basis functions from the cluster Hilbert space HIj, 
used in all possible combinations in (65). The motion of the clusters is 
described by a space Ma of functions of a system of relative coordinates ta 
between the cluster centers. 

Consider, for example, a 3-particle reaction XY + Z ^ X + YZ. Then 
^XY+z,o consists of products of approximate bound states il^xY and the 
ground state ipz] Mx+yz,o consists of products of the ground state ipx and 
approximate bound states ipyz, and M.xyz,q consists of sufficiently many 
states localized in the transition region to resolve the resonances of interest. 
Usually, one would keep the arrangements {XY, Z) and (X, YZ) in the re- 
duced description, and use the states belonging to the arrangement (XYZ) 
as doorway states for the transition regime. 
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Let C be the set of arrangements considered relevant for the reduced de- 
scription. The reduced multichannel state space is then the dual of the 

space 

Heff = 

or the properly symmetrized subspace in case of indistinguishable clusters. 
The inner product in is given by 

'y^eflf'0eflf= ^ / drA^AirAfi^AirA)- 
The embedding map P : E[*g H* is given by 

Pi^eS = ^ Pa'iI'A, 

where the kth component of Pa '■ (H^-*)* H* maps a function of ta to the 
kth basis vector of modified to have the corresponding dependence on 
the center of mass r^(x) of the cluster coordinates x, 

{PAll^A)k{x) = 1pA{rA{x))'^Ak{x). 

The transpose : H* ^ H*jf is given by 

{P^^)a = P^ip for all A e Ces- 
By construction, P^LP is a direct sum of contributions of the form 

P^LPa = La- Ha, 

where Ha is the free Hamiltonian for the motion of the cluster centers and La 

is a symmetric ua x n^-matrix. In a basis of cluster eigenstates with energies 
Ei, La is the diagonal matrix formed by the energy differences E — E^. 

The construction of (HqjPq) is completely analogous, using a larger set of 
arrangements and/or channels that contain the doorway states. 

Thus the calculations have the same complexity as those for the coupled re- 
action channel equations (or resonating group method), as described, e.g., by 
WiLDERMUTH & Tang [28]). However, the present scheme is more flexible 
in that it can incorporate information from doorway states. In principle, by 
increasing the size of the doorway state space, it is capable of arbitrarily 
accurate approximations to the full dynamics, and shares this feature with 
the two Hilbert space method of Chandler & Gibson [3, 4] and with a 
technique by Goldflam & KowALSKi [11]. 

To solve the reduced Schrodinger equation (and, if a similar construction 

is used for the doorway operator, the equations for the form factor), the 
whole arsenal of methods developed in the applications is available. Binary 
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arrangements can be handled by Lippmann-Schwinger equations (see, e.g., 
WiLDERMUTH & Tang [28], Adhikari & KowALSKi [1, Chapter 3]), and 
3-cluster arrangements by the Faddeev [6] connected kernel approach (see, 
e.g., Glockle [13, Chapter 3], KuKULiN et al. [17], Adhikari & Kowal- 
SKI [1, Chapter 7]). 
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